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Abstract 



We provide a new way of constraining the relative scintillation efficiency £ e g for liquid xenon. Using a simple estimate 
' | for the electronic and nuclear stopping powers together with an analysis of recombination processes we predict both the 
ionization and the scintillation yields. Using presently available data for the ionization yield, we can use the correlation 
^vq between these two quantities to constrain C e g from below. Moreover, we argue that more reliable data on the ionization 
i yield would allow to verify our assumptions on the atomic cross sections and to predict the value of £ e ff ■ We conclude that 
^ the relative scintillation efficiency should not decrease at low nuclear recoil energies, which has important consequences 
for the robustness of exclusion limits for low WIMP masses in liquid xenon Dark Matter searches. 



in 



Keywords: Dark Matter search, liquid xenon, WIMP, relative scintillation efficiency 



i-H 1. Introduction 
i 

O 
H 



CO 

c3 



(N 
> 
O 
OS 
OS 

m 



X 

S-H 

C3 



The recent results of XENON100 HJ |2] correspond to a 
major increase of sensitivity for the Dark Matter (DM) 
searches and further improvements are expected soon 
when more data will be released. However, the analysis of 
the detector properties turns out to be subtle for low recoil 
energies, corresponding to low mass DM particles (below 
lOGeV). The interest in this parameter region is addi- 
tionally heated by claims of DM observation j3J U] . The 
problem is the proper reconstruction of the nuclear recoil 
energy from the primary scintillation signal (SI) in liquid 
xenon in the limit of low recoil energies. Direct experimen- 
tal calibration is rather difficult for low nuclear recoil en- 
ergy and is prone to large systematic uncertainties, which 
led to mutually contradicting experimental measurements 



below 10 keV [5HT5] (see figure 5b I 



A theoretical treatment of the problem requires the de- 
termination of the scintillation yield for a slow moving 
xenon atom. The task can be roughly divided into three 
parts: the problem of ionization and excitation probabili- 
ties in the individual collisions of xenon atoms, the prob- 
lem of simulating the propagation through the media, and 
finally the problem of possible recombination of the pro- 
duced free electrons and ions. The outcome of such a theo- 
retical treatment would be predictions for both the scintil- 
lation and the ionization signals produced in the detector. 
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At sufficiently high energies the total electronic excita- 
tions in the atomic collisions can be reasonably well de- 
scribed by Lindhard's theory [16]. It approximates the 
process by point-like interactions between the incoming 
nucleus and electrons in the electron cloud of the tar- 
get, and is applicable for the case of nuclear recoil ve- 
locity v m 1=3 vo = e 2 /h (or E m f=s lMeV). However, for 
smaller energies, the individual collisions are much harder 
to calculate and require a nonperturbative analysis of the 
electronic movement |17j . For this reason, it is very diffi- 
cult to describe inelastic collisions at energies much below 
100 keV. 

In this article we do not attempt to make an ab-initio 
theoretical calculation of all the above processes, which is 
a very difficult task. Instead we make use of theoretical 
connections between the scintillation and ionization yields 
and the fact that the ionization yield is measured more re- 
liably at low nuclear recoil energies than the scintillation 



yield [T^] (see figure 5a I . For this purpose we will consider 
various possible modifications of the electronic stopping 
powers of liquid xenon using a simple parameterization in 
order to compare the resulting scintillation and ionization 
yields with experimental data. From the data for the ion- 
ization yield we observe that we cannot choose stopping 
powers that lead to small electron excitations at small re- 
coil energies without spoiling the fit of the ionization yield 
in figure |5a| This observation leads to the prediction that 
the scintillation yield also cannot decrease much at small 
nuclear recoil energies. 

The paper is organized as follows: in section [2] we review 
the process of generation of scintillation light and ioniza- 
tion in liquid xenon. Section [3] introduces the notion of 
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the nuclear and electronic stopping powers, which corre- 
spond to the analysis of individual xenon scattering events. 
Moreover, we calculate the total energy in the electronic 
excitations, and then analyze the recombination process 
in section U] In section [5] all the results are collected and 
translated into the ionization and scintillation yields and 
compared with the experimental data. Section [6] contains 
our conclusions and further prospects. 



2. Production of scintillation light in 
liquid xenon 

A large variety of effects must be taken into account to 
describe all physical processes that lead from the initial re- 
coil to the production of scintillation light in liquid xenon. 
Specifically, we expect the following steps [18_: 

• In a WIMP-like scattering event, an energy of 1- 
100 keV is transferred to the nucleus^ As the cor- 
responding recoil velocity is well below the Fermi ve- 
locity of the most loosely bound electrons, we expect 
the atom to remain neutral in the scattering process. 

• The recoiling atom will scatter off neighboring nu- 
clei. While most scattering events are expected to be 
elastic, there will occasionally be inelastic collisions 
leading to excitation or ionization of either (or both) 
of the atoms. 

• After each scattering process both atoms will continue 
their propagation with a fraction of the initial recoil 
energy. Consequently, both can again scatter elasti- 
cally or inelastically off other atoms. 

• During the process of thcrmalization the recoiling 
xenon atoms will leave behind a large number of ion- 
ized or excited xenon atoms — distributed along many 
branches of the initial track. 

• The free electrons will now either recombine with sur- 
rounding ions to form excited xenon atoms or escape 
from recombination. The fraction of escaping elec- 
trons will depend on the strength of the applied elec- 
tric drift field, but some electrons will escape even in 
the absence of a field. 

• Excited xenon atoms are free initially, but will soon 
be self-trapped and form excimers. These excimers 
emit vuv scintillation light on the transition to the 



2 Note that in the context of liquid noble gas detectors nuclear 
recoil energies are often quoted in keVnr. This unit is chosen to 
emphasize that any energy reconstructed from an SI signal depends 
on the effective scintillation yield (and is therefore not necessarily 
physical). Since we are concerned with actual physical processes and 
not the detector signals in this paper, we will give nuclear recoil 
energies in keV. 



ground state. In a simplified picture, the process is 
Xe* + Xe -> Xe? 



x e : 



2Xe + hi 



(1) 
(2) 



• In some cases, especially at high excitation density, 
two exited xenon atoms will combine to produce only 
one scintillation photon. This process, known as biex- 
citonic quenching, will effectively reduce the scintilla- 
tion yield. 



At first sight, the large number of steps seems to make it 
very hard to disentangle possible ambiguities. A decrease 
of scintillation efficiency at low recoil energies could, for 
example, be equally attributed to a decreasing cross sec- 
tion for inelastic scattering (for example due to threshold 
effects), a different track structure, an increasing fraction 
of escaping electrons or a stronger quenching mechanism. 
This ambiguity can be lifted at least partially by consid- 
ering not only the effective scintillation yield of nuclear 
recoils, but also the ionization yield (see also [H]). This 
quantity is much better measured but has been — to the 
best of our knowledge — ignored in all previous attempts 
to give a theoretical model for the scattering process in 
liquid xenon. 

The sum of ionization and scintillation, which we will 
refer to as the total electronic excitation, should corre- 
spond to the total energy lost in inelastic collisions. Con- 
sequently, it should only depend on the scattering cross 
sections and not on the processes occurring later, such as 
recombination, which will only lead to a redistribution be- 
tween ionization and scintillation. Thus, if both signals 
showed a similar energy dependence, this would suggest 
a general suppression of inelastic scattering at low ener- 
gies — which is what one might naively expect. However, 
one actually observes experimentally a strong increase of 
the ionization yield at low energies. This observation in- 
dicates that the suppression of the scintillation signal at 
low energies does not result from the actual inelastic scat- 
tering processes, but from the large number of escaping 
electrons. 

The effective scintillation yield of nuclear recoils is usu- 
ally described by the dimensionless quantity £ ff, called 
relative scintillation efficiency. It relates the SI scintilla- 
tion signal to the physical recoil energy of the nucleus E nr 
as 

Enr = 4 ■ (3) 

L y is the light yield for 122 keV electron recoils and S e ,n 
are the electric field quenching factors for electronic and 
nuclear recoils. Thus, C c g quantifies the suppression of 
scintillation for nuclear recoils compared to 122 keV elec- 
tron recoils at zero electric field. 
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3. Stopping powers of liquid xenon 

In this section we will describe the interactions of neutral 
xenon atoms and discuss possible scattering processes at 
energies of a few keV. The quantities we are interested in 
are the rate at which energy is transferred to recoiling nu- 
clei by elastic collisions and the rate at which electrons 
are excited by inelastic collisions. These quantities are of- 
ten called nuclear stopping power and electronic stopping 
power. 

3.1. Electronic stopping power 

The electronic stopping power is defined as the average 
energy which an atom loses to electronic excitations per 
distance travelled through the detector, (dE/ dx) . In a 
"semiclassical" approach, an electron is excited when it col- 
lides with a nucleus. The stopping power should therefore 
be proportional to the electron mass density no, their ve- 
locity v-p, the momentum transfer cross section cr tI (v-p), 
and the velocity of the incoming particle v. In fact, elec- 
tronic stopping is often described by 



(4) 



dE\ 

— J = n vv F a tr (v F ) 



We should take a moment to discuss the validity of such 
a semiclassical approach. Equation Q is in fact based on 
several assumptions: 

• Instead of describing the atomic system by a many- 
particle wavefunction, we claim that only the electron 
density is relevant to the problem, and we ignore the 
modification of the electron density during the atomic 
collision. 

• Collisions between electrons and nuclei are treated like 
classical point like interactions. 

• Electrons are assumed to be free. Consequently, no 
minimal energy transfer is required for an excitation. 

For a uniform electron gas the electronic stopping power 
is proportional to the velocity of the incoming particle, a 
result that has been derived in 1161. The authors obtain 



^) =V8ne 2 a ( ZN-- 
dx J v 



(5) 



where clq and vq are Bohr radius and Bohr velocity re- 
spectively and £o is an empirical parameter, ofterj^Jset to 
Co = Z 1//6 . N = 13.76 nm -3 is the number density of 
xenon atoms. This formula should also be valid for nu- 
clear velocities v < Vq. 

However, there is no reason that we can extrapolate the 
velocity proportional behavior all the way down to v = 0. 



In fact, at very low velocities, departures from velocity- 
proportionality have been observed experimentally [21 . 
Moreover, there are several theoretical arguments in fa- 
vor of a more rapid drop of ( dE/ dx) e for v < vq. The two 
most important ones are threshold effects and Coulomb 
effects. 

The argument for threshold effects essentially goes as 
follows: in an elastic collision between a nucleus and a free 
electron, only a fraction j m e /m^ of the nucleus energy can 
be transferred to the electron. For nuclei with energies in 
the keV range, the resulting electron energy is at most a 
few eV — so we can no longer ignore gap energies or the 
work function of xenon. There has been a long and intense 
discussion on whether such threshold effects are present 
or not (see for example [26 ). However, many experiments 
and theoretical considerations report electronic excitations 
far below the naive threshold |17l |2"7] . 

An effect from Coulomb repulsion is expected, because 
at very low relative velocities colliding nuclei will not pen- 
etrate the electron clouds of each other strongly. Conse- 
quently, with decreasing energy the recoiling nucleus will 
probe only regions of lower electron density [28] . 

Both of these arguments are doubtful, because they con- 
tinue to exploit the point-like interaction of the nucleus 
with the electron. However, for low nuclear velocity the 
electron clouds rearrange during the collision (or, in a more 
semiclassical language, the electron makes several rota- 
tions in the combined electric field of the two colliding 
atoms during the collision). This effect leads to a much 
more complicated non-perturbative mechanism of the en- 
ergy transfer to the electron. For an analysis of such col- 
lisions in the case of simple atoms, we refer to |17j . 

Ideally, exact or approximate quantum mechanical cal- 
culations for the Xe-Xe scattering process must be per- 
formed. Lacking such calculations, we will continue to use 
equation ([5| even for low energies. However, we will in- 
troduce a correction factor F(v/vq) to parameterize our 
ignorance of the cross section: 



— ^F(v/v ) — 



(6) 



In the end, we will compare our results with the exper- 
imental data on ionization and scintillation for different 
choices of F(v/vo). Until then, we will use F(v/vq) = 1 
unless explicitly stated otherwise. 

To conclude this section, we provide dimensionless quan- 
tities instead of (dE/ dx) e , which are usually preferred in 
the literature. Therefore, we define the reduced energy e, 
reduced distance p and a dimensionless electronic stopping 
power s e by 



c 
P 



E 



3 In an independent derivation, Firsov |25| obtained the same for- 
mula with £3 1.63. 
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where a = 0.626a Z~ 1 / 3 is the Thomas-Fermi screening 
length. For liquid xenon, the reduced energy e, which we 
will use throughout the rest of the paper, can be expressed 
as e = 1.05 • lCT 3 £ nr / keV. 

3.2. Nuclear stopping power 

The second quantity needed to calculate the amount of en- 
ergy lost to electronic excitations is the nuclear stopping 
power, corresponding to the probability for elastic scatter- 
ing of two xenon atoms. To calculate the cross section, we 
approximate the electron wave functions by the electron 
density, and ignore modifications of the electron clouds 
during the collision. The energy transfer in such a colli- 
sion depends in general on the energy E nr of the projectile 
and the scattering angle 0. However, it turns out that due 
to scaling properties all relevant functions depend only on 
the combined variable 129 



V 



eSm 2' 



(10) 



where e again denotes the reduced energy. 

The differential cross section for elastic scattering can 
then be written as 



da 
dr) 



(11) 



where a is again the screening radius. The function f(rj) 
depends on the screening function that we adopt to de- 
scribe the charge density. For a large number of screening 
functions, f(rj) can approximately be written as [30 
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From f(rj) we can obtain the dimensionless nuclear stop- 
ping power s n (e): 



Sn(e) = \ J dr//(r?) . 



(13) 



Both quantities, f{rf) and s n (e) have been calculated 
by various authors with differing results. These differ- 
ences, especially in the low energy region, are due to differ- 
ent approximations for the screening function. Lindhard 
et al. favor the Thomas-Fermi screening function corre- 
sponding to in — 0.333, q = 0.667, A = 1.309. However, 
today it is generally agreed that the Thomas-Fermi screen- 
ing function overestimates the potential at large distances 
and therefore gives too large stopping powers at low en- 
ergies [21]. One therefore often prefers the Moliere or the 
Lcnz- Jensen screening functions that show better agree- 
ment with experimental data. They correspond to the 
parameter choices m — 0.216, q = 0.570, A = 2.37 and 
m = 0.191, q — 0.512, A = 2.92, respectively. For a com- 
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(b) Nuclear stopping power 

Figure 1: Comparison of different choices for the nuclear 
stopping power. 

parison of the different screening functions, see figure [T| 

The most reliable way to calculate the screening func- 
tion is to use Hartree-Fock methods instead of an analyt- 
ical approximation. Ziegler et al. [3T] find the best agree- 
ment with experimental data for the (so-called universal) 
screening function given by the following expression]^] 



Sn(ez) 



ln(l + 1.1383e z ) 



2[ez 



0.013214 21226 



0.19593e£ 5 ] 



(14) 



From this, f(rf) can be calculated using f(x) = -xr [xs(a;)]. 

In contrast to the universal screening function from 
Ziegler et al. the advantage of the simpler screening func- 
tions is that f{rf) can be approximated by a power law for 
very small values of rj (i.e. 77 < 10~ 4 ): 



m - w 



(15) 



As we will not need this property, we will use the univer- 
sal screening function in the following, because it is the 



4 It is important to notice that Ziegler et al. use a slightly different 
definition for the reduced energy, because they assume a different 
screening length. For xenon, the conversion factor is ez = 1.068e. 
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most accurate at low velocities and agrees best with ex- 
perimental data. Note that the uncertainty related to the 
nuclear stopping power can be absorbed into the correction 
factor introduced in equation Q6J) , so that our discussion 
remains general. In figure [6j we will show how our final 
results would be affected by making a different choice for 
the screening function. 

3.3. Total electronic excitation 

In this section we will combine the nuclear stopping power 



s n (e) from equation ( 14) and the electronic stopping power 
s e (e) from equation (]4j) in order to predict what amount 
of the initial recoil energy is transferred to electronic ex- 
citations. We will denote the total energy in electronic 
excitations by /c(e) and also define the quotient 



£(e) = 



s„(e) 



(16) 



First of all we should try to discuss the general trends 
which we expect for k(c) and £(e). At energies above one 
MeV, inelastic collisions will dominate because the elec- 
tronic stopping power grows proportional to \/E, while the 
nuclear stopping power decreases after reaching its max- 
imum around 100 keV. Consequently, «(e) « e in this 
energy region. At energies of a few keV and below, the 
nuclear stopping power is much larger than the electronic 
stopping power. Consequently, we expect re(e) -C e and 
£(e)«l- 

In the low energy region, a first estimate for re(e) is given 

by 

e • s e (e) 



«(e) 



(17) 



s n (e) + s e (e) 

Making this estimate, we have included only electronic ex- 
citations produced by the primary recoiling nucleus. We 
expect, however, that at least some recoiling nuclei from 
secondary elastic collisions still have enough energy to in- 
elastically excite other atoms. Thus, part of the energy 
lost in elastic collisions can still be transferred into elec- 
tronic excitations. 

One way to take this effect into accout would be to per- 
form numerical simulations. Here, we will instead rely on 
Lindhard et al. |32j . who have derived an integral equation 
to determine n(e). Under the approximation that most 
electronic excitations occur at large impact parameter and 
have only small energy transfer, the authors show that for 
£(e) <C 1, k(c) oc e • £(e). Thus, we only need to modify 



our simple estimate in equation ( 17 1 by introducing a con 



stant of proportionality a in order to include the effect of 
secondary nuclear recoils: 



k(c) = ae£(e). 



(18) 



This constant of proportionality can be calculated ana- 
lytically, if £(e) can be described by a power law. As we 
want to consider cases where £(e) cannot be described by a 
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Figure 2: The function £(e) for different choices of the 
nuclear stopping power. The electronic stopping power 
was calculated using F(v/v ) = 1. 



power law, we cannot calculate a. However, we will show 
that this presently unknown constant can be absorbed into 
another quantity which we can determine from experimen- 
tal data. Since we expect the simple estimate in equa- 
tion (17 1 to underestimate k, we require a > 1, but still 
of order 1, for consistency. 



It should be emphasized that equation ( 18 1 does take 



into account secondary recoils produced along the entire 
track of the primary nucleus. Its simple structure is a 
result of the assumption e <C 1, which implies that the nu- 
clear stopping power is much greater than the electronic 
stopping power and both decrease rapidly as the recoiling 
nuclei loose their energy. More accurate ways to calcu- 
late k have been considered in the literature (see for ex- 
ample [33]). However, the original result from Lindhard's 
theory has the virtue that k oc F(v/vq), which will simplify 
the subsequent analysis. 

For Thomas-Fermi screening, we obtain from equa- 
tion (15 1 that £(e) oc e 017 and therefore n(e) oc e 117 . Con- 



sequently, Lindhard's theory predicts an increasing £(e) as 
the nuclear recoil energy increases from 1 keV to 100 keV. 
This result has often been quoted as a possible explanation 
for the energy dependence of the scintillation yield in liq- 
uid xenon. However, we understand now that this result 
strongly depends on our choice for the nuclear stopping 
power. As argued in section [3~2| choosing a Thomas-Fermi 
screening function will tend to overestimate the nuclear 
stopping power at low energies. Consequently, we must ex- 
pect to underestimate £(e). It appears much more reason- 
able to choose a nuclear stopping power that agrees better 
with experimental data (see also [33] )• In fact, choosing 
the universal stopping power from Ziegler et al. which is 
still closest to Thomas-Fermi, already changes the behav- 
ior of £(e) considerably (see figure [2]). Now, £(e) is no 
longer increasing monotonically, but develops a minimum 
around a few keV, remaining almost constant in most of 
the region we are interested in. 
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As we have seen, the uncertainties concerning the frac- 
tion of energy deposited in electronic excitations remain 
quite large. Although Lindhard's theory has been quoted 
frequently in the context of the effective scintillation yield, 
it appears difficult to obtain even a general tendency from 
this theory. We therefore believe that the underlying as- 
sumptions on the electronic and nuclear stopping powers 
are too weak to support the usual conclusion, that the 
energy dependence of £ e ff is due to a general suppression 
of electronic excitations at low energy. For this reason, a 
precise calculation of the relative scintillation efficiency is 
very difficult. To make progress, we will need to extract 
additional predictions from our theory that can be com- 
pared to experimental data to constrain our model. As 
we will show below, a reasonable agreement with experi- 
mental data is achieved for an almost constant £ (e) at low 
energies. 

4. Recombination 

Now that we have an estimate of the total energy in elec- 
tronic excitations, we need to determine how this energy is 
distributed between ionization and scintillation. This dis- 
tribution depends not only on the number of excited and 
ionized atoms produced initially, but especially on the re- 
combination rate. Recombination will occur whenever an 
electron and a ion produced in a nuclear recoil process ap- 
proach sufficiently close. The recombination rate should 
be proportional to the ionization density, which in turn 
is roughly proportional to the electronic stopping power 
s e (e). Consequently, we expect a higher ionization den- 
sity, and thus a higher recombination rate, at higher recoil 
energies. In this section, we will follow closely the Ph.D. 
thesis of Dahl [35] . 

As discussed in section [2] after all recoiling atoms have 
thermalized, we are left with a certain number of exci- 
tons, called N ex , and a certain number of ionized atoms, 
N{. We expect that a fraction r of the ionized atoms will 
recombine with free electrons, forming excitons that will 
eventually emit scintillation photons. The number of pho- 
tons produced should consequently be given bjj^j 



N 



P h 



N BK + r-N i = N i [r + 



Ni 



(19) 



We assume that the efficiency for the production of a 
scintillation photon in a recombination process is close to 
100%. Moreover, we assume that the fraction N ex /N is 
energy independent (see [3BJ f° r a discussion) , although it 
may depend on the nature of the recoiling particle^] This 



5 Strictly speaking, this equation describes the number of excited 
atoms after recombination. The actual number of photons produced 
will be reduced by quenching effects. However, for the moment, we 
neglect these effects, but will include them later on. 

6 For example, we expect N ex /Ni to be larger for the collision of 
two xenon atoms than for electron recoils, because the xenon atoms 
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Figure 3: Recombination fraction r as a function of the 
number of ionized atoms. Using the results from section [5] 
these values can be converted into reduced recoil energies, 
which are shown on the top of the graph. 



assumption is supported by computer simulations [35 (see 
also [40] for a discussion on ways to determine N ex /N ex- 
perimentally) Q Consequently, an energy dependence can 
only be introduced by the recombination fraction r. For 
later uses, we also define the number of electrons produced, 



N q = (1 - r) • 



(20) 



N cx and Ni are presently unknown. However, they 
should both be proportional to /c(e), which in turn was 
determined to be proportional to e£(e). As we kept the 
constant of proportionality undetermined, we can do the 
same thing for N cx + N , writing simply 



N cx + N l = Ni[l 



Ncx 

N 



= /?e£(e) 



(21) 



Here we have assumed that the mean energy for an exci- 
tation is approximately equal to the mean energy for an 
ionization, which means that N cx /N 1. We will see 
below, that this assumption is true to very good approxi- 
mation. Equation (21 1 allows to calculate N ex and Ni from 
k(c) once we have determined N ex /N and j3. If we also 
know the recombination fraction r(e), we can then infer 
N ph . 

The task is therefore to determine N ex /N and r(e) from 
an analytical model. In order to calculate the recombi- 
nation rate, one needs to describe diffusion processes for 
electrons and ions. In fact, various theories describing re- 
combination exist (for a review, see |42|V A model by 
Thomas and Imel [13] has been successfully used to de- 



can temporarily form molecular orbitals that enhance the probability 
for excitations I37H39I . 

7 An analysis of experimental data |41| . which appeared after this 
work was first submitted, further supports this assumption. 
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scribe recombination in liquid xenon |4"TI |4~4"] . This 
model gives 



10 



Ni 7^ V 4 



(22) 



where 7 is a free parameter of the theoryj^] 

The two free parameters 7 and N ox /Ni can be deter- 
mined either experimentally |H| or from Monte Carlo sim- 
ulations [85.. Here we take the values recently proposed 
by the XENON10 collaboration [H] as a conservative fit 
for the ionization yield in liquid xenon: 



N, 



1.09 , 



cx 

at 

7 « 0.032 . 



(23) 
(24) 



In |4T] it was argued that N cx /Ni has an uncertainty of 
about 15%. In addition. 7 can vary by as much as 20% for 
different values of the electric field [35] in such a way that 
the recombination decreases with increasing field strength. 
As these uncertainties do not significantly affect our con- 
clusions, we take the value from above at all electric fields 
in order to calculate the recombination fraction r accord- 
ing to equation ( 22 I (figure [3]) . 



5. Obtaining the scintillation yield 

Having calculated the total energy in electronic excitations 
and the recombination fraction, we are now able to predict 
both the ionization and the scintillation yield. However, 
one free parameter still remains in our theory: the propor- 



tionality factor f3 which we introduced in equation (21) 



Combining this equation with equation (22 I, we can write 



N q (e) 

me) 



imfi + ^p 
7 V 4 

1 + AW7V; ' 



(25) 
(26) 



iVq(e) has been measured experimentally (most recently 
in [T2l [T3]), so we can determine f3 from fitting equa- 
tion (25) to the available data. Instead of N q (e), one con- 



ventionally plots the ionization yield, which is defined as 
Q y {E nr ) — N q (E nr ) I E nl and measured in e~/keVnr. Al- 
though there is only one free parameter, a reasonable fit 
can be obtained setting (3 = (1.38 ± 0.10) • 10 5 (see figure 
[4J, which indicates that our description is sufficient 

In principle, a better fit could be obtained by taking 7 
as a free parameter and determine it from a fit to the data. 
In this case, larger values of both j3 and 7 are preferred, 
corresponding to a larger ionization yield at lower energies 



8 In fact, 7 can be expressed in terms of the recombination coef- 
ficient, a, the drift velocity, v, and the typical size of the track, a: 
7 = a/a 2 v. However, as these parameters are unknown, we may just 
as well take 7 as the parameter of the theory. 

9 The reason why /3 is so large is that fc(e) is a reduced energy. 



Manzur E = 1 kV/cm 
Theoretical Q„ 
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Energy / keV nl 



50 



100 



Figure 4: Experimental data for the ionization yield 
from |12| for an electric field of E = 1 kV/cm and the best 
fit of equation ([25} using N ex /Ni = 1.09 and 7 = 0.032. 



£(e) was calculated using F(v/vq) = 1 and the universal 
stopping power from Ziegler et al. 



and a steeper decrease towards higher energies. Doing so 
would also lead to larger values for £ e g at all energies, 
thus giving a rather optimistic prediction of the relative 
scintillation efficiency. We will therefore keep 7 = 0.032 
to give a conservative estimate of £ e ff, but point out that 
any attempt to improve the fit in figure [4] would actually 
lead to an increase of £ e g and not a decrease, as one might 
naively expect. 

We are now in the position to predict the relative scintil- 
lation efficiency. The value of N p ^ can be obtained directly 
from iV q because N q + jV p h = N cx + Ni = /3e£(e) and this 
sum is known once we have determined /3. In order to ob- 
tain £ e ff, we need to divide N p h by the number of photons 



-.rcf 
L ph ■ 



produced by the reference electron recoil at 122 keV, 
This value can be determined from the W p h(f3) value for 
xenon, which is the energy that an electron recoil must on 
average deposit in the detector to produce a scintillation 
photon^ 

Unfortunately, W p h(fl) is known only with some un- 
certainty. For electrons with an energy of 1 MeV, one 
observes W p h((3) — 21.6 eV [T^J. However, according 
to [45], the scintillation yield is energy dependent, with 
smaller values of T / F p h(/3) at lower energy of the electron 
recoil. At 122 keV, the scintillation yield is measured to 
be about 10% larger than at 1 MeV, corresponding to 
W p h(/?) ~ 19 eV. We will use this value for our analysis 
and therefore n p c ^ « 53ph/keV, but would like to empha- 
size that the uncertainty related to this value is quite large 



10 Note that Wph(/3) is different from W p h, which is the mean en- 
ergy required to produce either an excited or an ionized atom, cor- 
responding to the theoretical mean energy to produce a photon if 
there is full recombination and no quenching. W p h(/3) takes into ac- 
count these effects and is therefore larger than Wph- For a careful 
discussion of this issue see 1191 . 
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- at least 10%. 

Before giving our result for £ e ff we need to consider one 
more process, that has been neglected so far. As men- 
tioned in section [2] the number of excitons can be reduced 
by biexcitonic quenching (see [181 141)1 07]). The idea is 
that in collisions of two excited atoms, only one scintilla- 
tion photon is produced. Several authors have suggested 
a parameterization of this process in terms of Birk's satu- 
ration law |33| 148] . They introduce an energy dependent 
quenching factor q c \ given by 



1 



<7el = 



l + fc-s e (e) 



(27) 



where k is called Birk's constant and has been determined 
in [33J to be k = 2.015 ■ 10" 3 g/MeV cm 2 = 21.4 in re- 
duced units. The value suggested in [IS] is smaller by 
about 15%. Although it has been argued [UJ that the 



parameterization in equation ( 27 1 is not accurate in liq- 



uid xenon, the description is sufficient to note that biex- 
citonic quenching is only efficient at high recoil energies, 
when the density of excited atoms is large. From equa- 
tion ([9]), s e (e) = 0.166^/e, so k ■ s e (e) < 1 in the en- 
tire energy region considered and k ■ s e (e) < 0.25 below 
5keV. When only few excited atoms are produced, biex- 
citonic collisions are rare and cease to reduce the photon 



yield. Consequently, equation (27 1 should be a sufficiently 



good description in the energy region we are interested in. 
This statement is especially true if there is a mechanism 
that suppresses electronic excitations at low energy, corre- 
sponding to F(v/vq) < 1 for small v. 

Now we can write down our final result for C c g: 



..ref 



<Zcl(-Enr) 



(28) 



All parameters appearing in this equation have either been 
measured previously or were fixed above. Consequently, 
we can now plot £ e ff and compare it with available data. 
For F(v/vq) = 1 our model predicts a flat £ e ff at low recoil 
energies, giving roughly C c g = 0.09 at E nr — 2keV (see 
the red curve in figure 5b ) . Also, good agreement in both 
Q y and £ e ff is found between our results and the results 
that Sorensen obtained from Monte Carlo simulations of 
the nuclear recoil band |49j . 

We would like to emphasize that the Lindhard factor, 
meaning £(e), is still the dominating uncertainty of our 
model. In fact, we are under the impression, that this un- 
certainty has often been underestimated previously, since 
not even the general trend (suppression or enhancement) 
of electronic excitations at low recoil energies is fully clear. 
However, the point we would like to emphasize is that dif- 
ferent assumptions for the nuclear and electronic stopping 
powers affect the predictions for both ionization yield and 
relative scintillation efficiency in a correlated way. 

To illustrate this point, we also show in figure [5] ad- 
ditional curves corresponding to different assumptions on 




10 20 
Energy / keV nr 

(a) The ionization yield 



0.25 




0.15 



0.05 - 



10 20 
Energy / keV nr 

(b) The relative scintillation efficiency 

Figure 5: Predictions for the ionization yield Q y and the 
relative scintillation efficiency £ e ff compared to the ex- 
perimental data presented in |15| (and further references 
therein). The lines with different colors correspond to dif- 
ferent choices of the function F{v/vq): The red line, which 
we consider as the best description, corresponds to our pro- 
posal to modify Lindhard's theory by using the universal 
nuclear stopping power from Ziegler et al. (see also figure 
[2| and setting F(v/v ) = 1. The orange line has been 
obtained by assuming an enhancement of the electronic 
excitations at low energies as in equation ( 29 1 . The pur- 



ple line on the other hand corresponds to a suppression 
of the electronic stopping power by introducing a smooth 
cut-off for £ (e) as in equation ( 30 1 . While the enhance- 



ment is excluded by the data for the relative scintillation 
efficiency, in the case of a suppression the ionization yield 
becomes clearly inconsistent with data, thus limiting the 
scintillation efficiency from below. For all plots we have 
used P = 1.38- 10 5 . 
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£(e) (corresponding to different choices of the correction 
factor F(v/vq) in equation Q). In detail, the cases con- 
sidered are 

F cnh {v/v ) = l + exp(-50e) (29) 
F supp {v/v ) = i(l + tanh(50e- 0.25)) (30) 

corresponding to an enhancement and a suppression of the 
electronic excitation at low recoil energies, respectively^] 
F SU pp{x) is similar to the result obtained by Tilinin |20J in 
an attempt to include Coulomb effects in the calculation 
of the electronic stopping power. It could, in principle, 
account for the drop of the relative scintillation efficiency 
observed experimentally in [13J. -FcnhOz) is harder to moti- 
vate from a theoretical perspective, but could result from 
smaller nuclear stopping power at low energies. 

Due to recent measurements of £ e ff |14| , we can clearly 
exclude the possibility of an enhancement of the electronic 
excitations at low energies as described by F fm i 1 (x). How- 
ever, the remaining two predicted curves for C e s are sim- 
ilarly compatible with experimental data. From the the- 
oretical point of view, it would not be easy to justify a 
definite preference for one of these choices either. How- 
ever, choosing F(v/vq) < 1 will also significantly reduce 
the ionization yield. If £(e) decreases with decreasing re- 
coil energy (as it would be the case in the presence of 
threshold effects or Coulomb effects), one cannot account 
for the increasing ionization yield that we observe exper- 
imentally Our central observation is therefore that any 
model attempting to explain the scintillation yield of liq- 
uid xenon, must at the same time explain the ionization 
yield. Consequently, a general suppression of electronic ex- 
citations at low recoil energies is clearly incompatible with 
experiments. In contrast, a nearly constant value of £(e) 
at low energies can accommodate (and even predict) the 
opposing trends seen in scintillation and ionization yield, 
because of the energy dependent recombination fraction. 

There are several sources of errors in the present anal- 
ysis. The main theoretical errors come out of the calcu- 
lation of the nucelar and electronic stopping powers and 
the understanding of the ionization process. The first, cor- 
responding to the nuclear stopping power, is comparably 
small, as far as the knowledge of elastic scattering is quite 
good. To emphasize this point, we also show in figure [6] the 
effect of taking different nuclear stopping powers to calcu- 
late the electronic excitations £(e). Indeed, the universal 
stopping power from Ziegler et al. is not only best moti- 
vated theoretically, but also gives the best description of 
both ionization and scintillation. The uncertainties in the 
electronic stopping powers and ionisation process, on the 
other hand, can change the result considerably. Lacking 
a reliable calculation of xenon inelastic scattering at low 
energies, we cannot provide an exact estimate of the in- 
duced theoretical error — the wide spread of the graphs in 



11 Note that e is of course a function of v/vq. 
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Figure 6: Predictions for the ionization yield Q y and the 
relative scintillation efficiency C e g compared to the ex- 
perimental data presented in |15| (and further references 
therein). The lines with different colors correspond to dif- 
ferent choices of the nuclear stopping power as in figure [2j 
For all plots we have used (3 = 1.38 • 10 5 . 



Fig. [5] can be considered as a spread of theoretical predic- 
tions. However, after choosing the value for £(e) from the 
ionisation yield measurements (upper plot), the respective 
curve on the lower plot will have much smaller error. 

Our results rely on the assumption that the ratio N ex /N[ 
does not vary strongly with energy. If we allow for an ar- 
bitrary energy dependence of this quantity and simultane- 
ously vary £(e), we could obviously fit any measurement 
for the ionization yield and the relative scintillation effi- 
ciency. However, from all available data, such an energy 
dependence is not expected. Moreover, to suppress scin- 
tillation and enhance ionization at low energies, N cx /Ni 
would have to decrease, implying that excitation becomes 
less likely compared to ionization. Such a behavior would 
most likely contradict the fact that less energy is required 
to excite a xenon atom than to ionize it. 
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6. Conclusions 

We studied in this paper the behavior of low energy nuclear 
recoils in liquid xenon which is important for Dark Matter 
searches. An ab-initio analysis of the problem would re- 
quire a quantum mechanical calculation of exclusive cross 
sections in Xe-Xe collisions (elastic, ionizing and leading 
to excited atoms) and a Monte Carlo simulation of the 
propagation of the nuclear recoil and further recombina- 
tion. The collision analysis should be made for low impact 
velocities, taking into account the electronic structure of 
the outer shells of xenon (for example, using Density Func- 
tional Theory). Alternatively, these cross sections could 
be determined in experiments on scattering of individual 
xenon atoms. 

To do so is obviously a formidable task and we pursued 
therefore a simplified effective framework. Specifically we 
combined Lindhard's theory for the initial production of 
electron excitations in atomic collisions with the well moti- 
vated assumption of energy independent partition between 
ionization and excitation in the underlying collisions. In- 
cluding an analysis of the recombination processes we ob- 
tained both the ionization and scintillation yields which 
show a correlated functional behavior when the stopping 
powers at low energies are varied. 

We argued that this correlated behavior allows an inter- 
esting consistency check when ionization and scintillation 
are measured simultaneously. We also argued that this 
correlation allows to use low recoil energy data for ioniza- 
tion to predict the low recoil energy dependence of scintil- 
lation. Using existing low recoil energy data for ionization 
we showed that it favors a constant behavior of C c g for low 
energies and that it allows to exclude the possibility of C e f{ 
dropping rapidly to zero below some threshold of 5 keV or 
more. In the analysis of the first XENON100 data [T], 
a constant behavior of C e g at low energies was assumed, 
which led to discussions in the literature about the reliabil- 
ity of the bounds for low WIMP masses. Our study shows 
that a constant £ e g is strongly favored by the ionization 
data and this strengthens therefore the low WIMP mass 
bounds found by the XENON100 collaboration. 
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